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LONGITUDINAL STABILITY, CONTROL EFFECTIVENESS, 

AND DRAG CHARACTERISTICS AT TRANSONIC SPEEDS OF A ROCKET- 
PROPELLED MODEL OF AN AIRPLANE CONFIGURATION HAVING AN 
UNSWEPT TAPERED WING OF ASPECT RATIO 3.0 
AND NACA 65A004.5 AIRFOIL SECTIONS 


By John C. McFall, Jr. and James A. Hollinger 
SUMMARY 


A flight investigation has been conducted on an airplane configura- 
tion mounting an unswept tapered wing of aspect ratio 3.0 with NACA 65A004.5 
airfoil sections. Static and dynamic longitudinal stability, control effec- 
tiveness, and drag characteristics were obtained. Comparisons are pre- 
sented between the model in the present investigation and similar models 
mounting wings of the same plan-form dimensions but heaving hexagonal air- 
foil sections. 


For models having wings of the same material, the NACA 65A004.5 air- 
foil model had lower values of lift-curve siope than that of the hexagonal 
airfoil section model. Some nonlinearity in the static stability of the 
complete configuration was observed at high subsonic speeds. The sharp 
static stability variations for all models at Mach numbers. between 0.90 . 
and 0.95 were indicated to be plan-form effects. The portion of the total 
airplane normal force carried by the exposed wing was approximately con- 
stant over the Mach number range from 0.73 to 1.33. Some buffeting of the 
model was indicated at high subsonic speeds at lift coefficients only a 
little below maximum. Тһе zero-lift drag of the wing on this configura- 
tion was only a small part of the total model drag. The drag rise of the 
model occurred at a Mach number of approximately 0.9. 


INTRODUCTION 


The rocket-propelled airplane configuration model flown in this 
investigation is the third model, in a general longitudinal-stability 
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research program, mounting a wing of aspect ratio 3.0, with 16° swept. 
quarter chord and taper ratio 0.40. The wing in this investigation dif- 
fered from the two models previously flown by having NACA 65A004.5 air- 
foil sections rather than hexagonal airfoil sections (ref. 1). The data 
obtained were analyzed by the methods presented in the initial report of 
this program (ref. 2). The results presented are those obtained from 
analysis of the model response to deflections of the horizontal tail 
which was moved in an approximate square-wave program. The longitudinal 
stability, control effectiveness, and drag characteristics of the con- 
figuration were determined from the flight time histories for a Mach 
number range of 0.73 to 1.33. The model was flown at the Pilotless Air- 


craft Research Station, Wallops Island, Va. = 


SYMBOLS 
Cy normal-force coefficient (ва/ a) 
Сс chord-force coefficient CL “> 
Ст, lift coefficient (CN cos а - Сс sin a) 
CD - drag coefficient (Сс cos a + Cy sin a) 
Cmtotal total pitching-moment coefficient a | © 
Оу wing normal-force coefficient  - um Ap v 
CN complete-model normal-force coefficient 
an normal acceleration as obtained from accelerometer, | 
ft/sec? : 
а; longitudinal acceleration as obtained from.accelerometer, 
ft/sec B 
g acceleration of gravity, #%/весе 
i distance between accelerometers, ft 
р free-stream static pressure, lb/sq ft : Е 
+ distance between model c.g. and tail center of pressure, ft 
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Po standard sea-level static pressure (2116 1b/sq ft) 

с wing mean aerodynamic chord, ft 

be/2 exposed wing semispan, ft 

е wing chord, ft 

G torsional modulus of elasticity, lb/sg in. 

y lateral distance from fuselage center line, ft 

mg wing torsional stiffness parameter, in.-lb/radian 

m couple applied near wing tip, in.-lb 

Ө local wing twist angle produced by m, radians; or angle 


of pitch, deg 


a angle of attack, deg 

5 elevator deflection, deg 

M Mach number 

S wing area (including the area enclosed within the 
fuselage), sq ft | 

W weight, 108 

q free-stream dynamic pressure, lb/sg ft 

ly moment of inertia in pitch, slug-ft^ 

P period of oscillation, sec 

R Reynolds number, based on wing mean aerodynamic chord 

T1/2 time to damp to one-half amplitude, sec 

д pitching acceleration, deg/sec? 

y flight-peth angle, deg 

y . ay/at 
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subscripts: 

T trim | | 
5 (ао/аъ) (2/27) (1/57. 3) | : U 
q (аө/аъ) (c/2V) 

W wing = i 
A complete model 

F fuselage | 

5 tail i 


MODELS AND APPARATUS 


Models zu = 


The model used in this investigation is the third in a general 
research program to investigate thin unswept wings of aspect ratio 3. 
This model is referred to in this paper as model C. The general model | 
dimensions are presented in figure 1, and three photographic views of the 
configuration are shown in figure 2. For comparison purposes, information 
is included herein on the two thin unswept-wing models from reference 1, 
which are referred to as models A and B. The three models were closely 
similar and differed only in the following respects: wing airfoil and 
material, fuselage length, and vertical tail. 


Models A and B had wings of hexagonal section and model C had a wing 
with МАСА 65A004.5 airfoil sections. The material of construction of the. 
wings was steel for model A and aluminum for models B and C. Model C had 
two inches of additional fuselage length in the cylindrical part of the 
fuselage ahead of the wing. The tail of model U had an extra vertical 
fin added below the fuselage to prevent the erratic behavior of longitu- 
dinal oscillations mentionéd in the analysis section of reference 3. 

The all-movable horizontal tail was deflected from 1.10° to -1.95° in an 
approximate square wave pattern. A detailed description of models A 

and B is presented in reference 1. Model C weighed 131.2 pounds, had a 
moment of inertia about the axis of pitch of 12,84 slug-feet square, and 
had its center of gravity at 12.5 percent of the mean aerodynamic chord. 
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Instrumentation 


The model was equipped with a 10-channel telemeter to transmit data 
from the model during flight. The continuously recorded measurements 
were as follows: two normal accelerations, longitudinal acceleration, 
transverse acceleration, angle of attack, two wing normal-force measure- 
ments, control position, total-head pressure, and static pressure. The 
normal accelerations were measured at the center of gravity and at the 
nose to provide a measurement of pitching acceleration. The wing was 
mounted on a beam-type balance from which wing normal forces were measured. 
One wing normal-force gage was an inductance type of gage used previously 
on another model (ref. %) and the other was a strain-gage type installed 
for development purposes only. A vane-type angle-of-attack indicator was 
mounted on a sting protruding from the nose of the model (ref. 5). The 
total-pressure tube was located on a small strut below the fuselage. The 
static pressure was measured from a calibrated orifice located 4.9 inches 
behind the beginning of the cylindrical portion of the fuselage. 


Radar tracking units were used to obtain model range and elevation, 
and Doppler velocimeter radar recorded velocity against time. Radiosondes 
measured atmospheric conditions at the time of the flight. Fixed and 
manually operated 16-millimeter motion-picture cameras were used to pho- 
tograph the launching of the model and the first portion of the flight. 


Preflight Measurements 


The torsional-stiffness characteristics of the wing were determined 
by applying a couple at the tip and measuring the twist at several sec- 
tions. The results are shown in figure 3. The factor plotted, Gctip)/mg, 
is a nondimensional parameter which makes the result a dìrect function of 
wing shape only. For use in comparing the aeroelastic properties of this 
wing with other results the value of free-stream static pressure, as 
obtained in the model flight, divided by standard sea-level pressure is 
presented in figure № as a function of Mach number. 


For use in analyzing the test results for buffet information the 
model was subjected to vibrations of known frequency, in order to obtain 
the response of telemeter instruments in the model and vibrational char- 
acteristics of major components, such as wing, tail, or nose. 


Measurements were made of the weight of the moving parts of the wing 
balance and wing panels. These measurements were used in applying a cor- 
rection for inertia effects on the wing-balance readings by the method of 
reference В. 


At values of wing normal force above 800 pounds the wing root dis- 
placement was stopped by the fuselage because of fuselage flexibility. 
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This resulted in incorrect values of wing normal-force readings obtained 
from the deflection of the wing-balance inductance-gage instrument. For 
this reason, the values of wing normal force obtained from high-lift 
oscillations were not used for the wing-alone lift-curve slopes nor for 
the ratios of wing lift to airplsne.lift. Any reference to wing-alone 
data in this report includes interference. 


TESTS AND ANALYSIS 


Tests 


The model was leunched at an elevation of approximately 60° from a 
hydraulically elevated mobile launcher as shown in figure 2(c). The 
model contained no rocket motor and was boosted to maximum velocity by a 
6-inch-diameter ABL Deacon solid-fuel rocket motor. Drag separation of 
the model from the booster followed booster burnout by reason of the dif- 
ferent drag-weight ratios of the model and booster. 


The scale of the test is shown in figure 5(a), where Reynolds number 
is plotted as a function of Mach number. 


Analysis 


The data were analyzed for the decelerating portion of the flight 
after model separation from the booster. A small corréction for rate of 
pitch was applied to the indicated angles of attack to convert them to 
angles of attack at the center of gravity of the model (ref. 5). The 
wing normal-force measurements were corrected for inertia effects by the 
method of reference 1. 


Басһ abrupt movement of the horizontal tail produced a short-period 
oscillation in аһ, 0, and aj, which damped during the time that the 
horizontal tail remained fixed, and from these quantities Сі and Cp 
were reduced. Analysis of the period, damping, and trim angle of attack 
using the equations of motion for two degrees of freedom yields the 
desired aerodynamic parameters. A detailed presentation of this method 
of analysis is made in appendix A of reference 2. 


The additional instrumentation used in this investigation but not 
used in reference 1 was a wing normal-force balance and a normal accele- 
rometer in the nose. The values from the wing balance were used for 
obtaining wing-alone normal-force coefficients and wing contribution to 
the total lift. The normal accelerometer in the nose was used, together 
with the normal accelerometer at the center of gravity, to obtain Pe 
acceleration used in the following eguation: 
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where 
6 = 5 (а/е) nose - (an/&)cg] 


The data were recorded continuously on the telemeter record for all 
channels but, for the purpose of plotting and reduction of data, readings 
of the telemeter record were made at intervals of 0.01 second. А11 slopes 


were taken from plots of the first 1 or 15 cycles of each oscillation. 


Periods, damping, and trim values were determined from plots of entire 
oscillations. 


ACCURACY 


Previous experience has shown the accuracy of telemetered measure- 
ments to be 2 percent of the full-scale calibrated range of the instru- 
ment. The following table gives estimated values of the possible system- 
atic errors in the absolute values of Сул; Cc, and Суу, as affected by 


the instrument calibration ranges. 


1.30 6470 


1.10 4620 
.80 2240 


А consideration of all the factors involved indicates that the Mach 
numbers are accurate to +1 percent at supersonic speeds and +2 percent 
at subsonic speeds. Further errors in the aerodynamic coefficients may 
arise from possible dynamic-pressure inaccuracies which are approximately 
twice es great as the errors in Mach number. 


The errors in the measured angles of attack and horizontal-tail 
deflections, being independent of dynamic pressure, are not likely to 
vary with Mach number. The horizontal-tail deflections are estimated to 
be accurate to 30.19 and the increments in angle of attack to 30.29, 


—— 
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RESULTS AND DISCUSS ION 


Lift Limit and Buffeting 


The lift information and range of buffeting, analyzed by using 
vibrational characteristics of the model, are summarized in figure 5(b). 
Some additional buffet information at lower Mach numbers is presented in 
figure 5(c) as angle-of-attack variation rather than lift coefficient 
because of the low values of normal force and the consequent decreased 
accuracy of the lift coefficient at the lowér Mach numbers. A similar 
variation with Mach number of the boundary for unsteady air flow to that 
experienced by the model between М = 0.6 and М = 0.7 was found for 
thin airfoils in the study of fluctuating pressures on two-dimensional | 
airfoils in reference 6 апа was attributed to leading-edge flow attachment. 


Lift TM | = 


The lift curves obtained over the Mach number range from M = 0.7 
to M- 1.3 for this investigation are presented in figure 6. The data 
points shown are those reduced from the telemeter record at time intervals 
ОҒ 0.01 second. Тһе lift coefficients are plotted against angle of attack 
for the first 15 cycles of each oscillation. Difference in lift at a | 
given angle of attack in figure 6 is shown by symbols representing Е 
increasing and decreasing angie of attack. This difference was found to 
be of smaller magnitude for the round-nose airfoil in the present inves- 
tigation than for models A and В in reference 1. 


The average lift-curve slope of model С is shown in figure 7 plotted 
against Mach number along with the lift-curve slopes of models A and B 
(ref. 1). The data points shown in figure 7 represent points from the 
faired lines presented in figure 6. The variation of lift-curve slope 
with Mach number was similar for all models, with model C having a lower 
value than model B which was of the same material of construction. 


The lift effectiveness of the all-movable horizontal tail ЛОт/Лё, 
presented in figure 8 for model C, are in good agreement with the values 
given for models A and B, which were shown to be of the right order of 
magnitude in reference 1. The break in the curve for model C results 
from absence of trim data at the trim change for the high-lift tail 
setting. 


Wing normal-force coefficients plotted against total configuration 
normal-force coefficients during low-lift oscillation are presented іп _ 
figure 9. The high-lift. oscillation data are not presented for the rea- 


son given in the section entitled “Preflight Measurements." The wing 
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normal-force coefficients were determined from the telemetered data at 
O.Ol-second time intervals. Some difference is noted for increasing or 
decreasing angle of attack. 


The rate of change of wing normal-force coefficient with total con- 
figuration normal-force coefficient is presented in figure 10 plotted 
against Mach number. The values obtained indicate that the wing contri- 
bution to the total lift is approximately constant over the Mach number 
range M=0.7 to M= 1.3. The ratio of exposed to total wing area 
compared with the values of абу / 4 Ср indicates that the fuselage lift 


: is in proportion to the wing area enclosed within the fuselage. The 


value dCNg/dCNyp was obtained by using the values of асқу/4Сңд from 
figure 9 where dCy, has been corrected for tail lift (ref. 1) which is 
shown as асц,/а Си - A comparable vaiue of асқу/асқұт may be found іп 
reference ( where for a wing of aspect ratio 3 a similar wing-body com- 


CN 
bination at a Mach number of 0.25 had a value of Gy SW = O, Th. 
cawr 


Drag 


Figure ll presents the minimum drag coefficients of model C, and 
the round-nose and hexagonal-airfoil section drag models (ref. 8). The 
drag rise for model С occurred at a Mach number of 0.9. The drag for 
models A and B is not presented, since the configurations differed from 
Model C by the absence of the additional vertical fin and by having dif- 
ferent wing airfoil sections. A comparison of model C with the wingless 
model in reference 1 shows that the wing contribution to the total drag 
is only a small part of the total drag of the configuration. The differ- 
ence in the minimm drag coefficients for the drag models (fig. 11) indi- 
cates an appreciable effect of airfoil shape. A comparison of the effect 
of lift on drag (fig. 12) shows lower values for model B when compared 
with model C as might be expected from the higher values of lift-curve 
slope for model B. No difference in chord force with angle of attack is 


indicated by the values of асг/аст2, since the values are nearly equal 
1 
to the value for no leading-edge suction, = ===>, 
© се | 21«3 Ста, 
The maximum lift-drag ratios for models B and C are plotted against 
Mach number in figure 13. Also presented in figure 13 are the lift-drag 


ratios obtained by using minimum drag values including induced drag from 
the drag models (ref. 8). The magnitude of the (1/0) ах values obtained 


by using the minimum drag of the drag models (ref. 8) is representative 
of an airplane with a very low fuselage drag. 
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In figure 14 are presented the lift coefficients at which maximum 
lift-drag ratios occur for models B and C. | 


Static Stability 


The measured periods of oscillations of angle of attack for model C 
is shown in figure 15. The data converted to the static stability derive- 
tive Cy, at a center of gravity of 12.5 percent of the mean aerodynamic 


chord and (Cm, for models A and B at the same center of gravity are pre- 
sented in figure 16. A similar variation of “Cy, with Mach number of 


models A, В, and С indicates that the sharp variation between М = 0.90 
and М = 0.95 were plan-form effects. Above М = 1.0 model C had a 
smoother variation of Cm, with Mach number than models A and B. 

The variation of the aerodynamic-center location with Mach number 
for models A, B, end C is shown in figure 17. The same trend is observed 
for all three models with model C having a smoother variation with Mech 
number above М = 1.0 when compared with models A and B. | 

The values of pitching-moment coefficient obtained by using the two 
normal accelerometers presented in figure 18, have been corrected for 
damping by the following formula: 


Cm(9) = Cmtotai - (Сид + Спа) & - Сша? 
where 


i 


and Суб + Cm, У, 9, Спњођај and Сш are obtained from the flight 


data. The faired lines shown are those obtained from the readings of the 
periods. Some of the pitching-moment-coefficient curves when plotted 
against time did not oscillate about a mean value of zero; probably 
because of small errors in the readings of either or both of the normal 
accelerometers. When this situation occurred, the curves were shifted 

to produce a mean value of zero, as required by the equilibrium conditions. 


Damping in Pitch . zs 4 | БЕ 
The time to damp to one-half amplitude plotted against Mach number 


is presented in figure 19. The data converted to damping coefficient 


Ono + Cms, and plotted against Mach number are shown in figure 20 for 


о“ “УУ 
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models А, В, апа С. No difference in damping coefficient for different 
tail deflection was noted for models A and B. The reason for the disa- 
greement in damping coefficient between models A and B and model C is 

not known but the same sharp variation st М = 0.8 to М = 0.9 has been 
observed in similar configurations with different wings (refs. 3 and 4). 
These differences between models were not observed in static stability. 


Longitudinal Trim and Control Effectiveness 


The variation of trim lift coefficient and the trim angle of attack 
with Mach number for the two control positions are shown in figures 21(а) 
and 21(b), respectively. No abrupt changes were apparent in either 
quantity with Mach number. 


The variation of the change in trim angle of attack per degree 
change in horizontal-tail deflection with Mach number is presented in 
figure 22 for models A, B, and C. Good agreement was obtained among 
the models for this quantity as well as for change in trim lift coeffi- 
cient per degree change in horizontal-tail-deflection variation with Mach 
number shown in figure 23. 


CONCLUSIONS 


The flight investigation over the Mach number range from 0.75 to 1.3 
of an airplane configuration having a thin unswept wing of aspect ratio 3 
and NACA 65A004.5 airfoil sections indicated the following conclusions: 


1. The variation of lift-curve slope was similar for the model having 
an NACA 65A004.5 airfoil section and similar models with hexagonal airfoil 
sections. The model with the NACA 65A004.5 airfoil section had a lower 
value of lift-curve slope than a model with a wing of the same material 
but with hexagonal airfoil section. 


2. The part of the total airplane normal force carried by the 
exposed wing was approximately constant over the Mach number range from 
0.73 to 1.23. The fuselage lift was in proportion to the wing area 
enclosed within the fuselage. 


3. Buffeting of the model occurred at high subsonic speeds at lift 
coefficients only a little below the maximun. 


4. The zero-lift drag of the wing of this configuration was only a 


51811 part of the total-model drag. The drag rise for the model occurre 
at a Mach number of approximately M = 0.9. a 


aM 
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5. The static stability of the complete configuration varied with 
lift coefficient at subsonic speeds. A similar variation of static sta- 
bility with Mach number for both hexagonal and NACA 65A004.5 airfoil 
section wing models indicated that the sharp variations at Mach numbers 


between 0.90 and 0.95 were plan-form effects.— 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, | ee 


Langley Field, Va, 
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Figure l.- General arrangement of model. А11 dimensions in inches. 
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(b) Top view. 
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Figure 2.- Photograph of configuration investigated. 
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(с) Model on launcher. 


Figure 2.- Concluded. 
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Figure 3.- Measured torsional rigidity of the wing. 
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Figure 4.- Variation of static-pressure ratio with Mach number. 
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(a) Reynolds number based on wing mean aerodynamic chord. 
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Figure 5.- Range of tests. 
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(c) Variation with Mach number of angle of attack at which buffeting 


12 


PTET (11411 
_ LLL LEAL LLL UL 
Es И 


72089 22229 265 

92 

ССС 
6 7 8 9 Ke 


M 


first occurred. 
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(b) в = 1.10. 


Figure 6.- Variation of lift coefficient with angle of attack. 
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Figure 8.- Lift effectiveness of horizontal tail. 
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Figure 9.- Variation of wing normal-force coefficient with airplane 
normal-force coefficient during Several typical oscillations. 
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Figure 10.- Rate of change of wing normal-force coefficient with airplane 
normal-force coefficient. 
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Drag models( Reference 8) 
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NACA 65A004.5 airfoil 
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Figure 12.- Effect of lift on drag. 
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Figure 14.- Lift coefficients for maximum lift-drag ratio. 
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Figure 16.- Static stability of complete models. Center of gravity at 
12.5 percent mean aerodynamic chord. 
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Figure 17.- Aerodynamic-center location. 
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(b) 8 = -1.95. 


Figure 18.- Variation of pitehing-moment coefficient with lift coefficient. 
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Figure 20.- Damping derivative. 
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(b) Trim angle of attack. 


Figure 21.- Trim characteristics. 
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Figure 22.- Change in trim angle of attack per degree change in horizontal 
tail deflection. 
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Figure 23.- Change in trim lift coefficient per degree change in horizontal 
tail deflection. 
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